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We shall consider domain walls in a relativistic field-theoretical model 
defined by the following Lagrangian 

where $ is a single real scalar field, (?y M1/ )=diag(-l, 1,1,1) is the metric in 
Minkowski space-time, and A, M are positive parameters. 

Euler-Lagrange equation corresponding to Lagrangian (1) has the partic- 
ular time-independent solution 

3 

$(^)=$ tanh(^) (2) 

which describes a static, planar domain wall stretched along the x 3 = plane. 
Here $o = M/ (2\/A) denotes one of the two vacuum values of the field $ - 
the other one is equal to — $o- The parameter M can be identified with the 
mass of the scalar particle related to the field $. The corresponding Compton 
length l = M~ x gives the physical length scale in the model. Energy density 
for the planar domain wall (2) is exponentially localised in a vicinity of the 
x 3 = plane. The transverse width of the domain wall is of the order 2Zo- 

The issue is time evolution of a non-planar domain wall. Such domain 
walls can be infinite, consider, e.g., locally deformed planar domain wall or 
a cylindrical domain wall. They can also be finite closed, e.g., like a sphere 
or a torus. We shall restrict our considerations to a single large, smooth 
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domain wall. Such a domain wall is denned by a set of conditions which 
provide Lorentz and reparametrisation invariant formulation of the heuristic 
requirement that IqRi^ <C 1 at each point of the domain wall, where i?i,i?2 
denote local main curvature radia in a local rest frame of the considered 
infinitesimal piece of the domain wall. In this case it is possible to develop 
the presented below perturbative approach to the dynamics of the domain 
wall. In other cases, e.g. when Iq/Ri ~ 1 at certain points of the domain 
wall, the only practical tool is numerical analysis. 

Physical idea underlying the expansion in the width is rather simple. The 
domain wall is a stable solitonic object. One expects that each piece of the 
large, smooth domain wall in the local co-moving reference frame does not 
differ much from the planar domain wall. The differences, which are due to 
curvature of the domain wall, can be calculated perturbatively. 

The tricky point is to find a convenient reference frame co-moving with 
the domain wall. From the seminal papers [[I], |2|], which were devoted to 
dynamics of vortices but this does not make an essential difference here, we 
know that the co-moving reference frame should be based on a relativistic 
membrane. In the leading approximation [0], |2[ the membrane is of Nambu- 
Goto type, and at each instant of time it coincides with the surface on which 
the scalar field <3> vanishes. In the following we will call this surface the core 
of the domain wall. For instance, in the case of the static domain wall (2) the 
core is given by the plane x 3 = 0. When calculating corrections to the leading 
approximation one can still adhere to this identification of the membrane with 
the core, see, e.g., papers || [|, [J. The price for this is that the membrane 
is no longer of the Nambu-Goto type. Effective action for such membrane 
contains higher derivatives with respect to time, and probably is nonlocal |J. 
The reason for these unpleasant features is that the points at which the scalar 
field vanishes do not constitute a physical object (the physical object is the 
domain wall itself and not the core), and therefore the core, being a purely 
mathematical construct, can have strange from the physical point of view 
equation of motion. In papers |7j § we have proposed to use a membrane 
which is of Nambu-Goto type to all orders. This membrane coincides with 
the core only at the initial instant of time. 

Our description of the dynamics of the domain wall |7|, || involves the 
Nambu-Goto membrane and certain additional 2+1-dimensional fields de- 
fined on the membrane. They obey nonlinear equations of motion which are 
second order partial differential equations, and they describe time evolution 
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of the domain wall completely. Transverse profile of the domain wall, that 
is dependence of the scalar field $ on a variable changing in the direction 
perpendicular to the domain wall, is uniquely and explicitly given once the 
evolution of the Nambu-Goto membrane and of the 2+1-dimensional fields 
is found. The dependence on the transverse variable is calculated perturba- 
tively, and one should observe that it is the case of a singular perturbation. 
Our approach to constructing the perturbative expansion follows a method 
used in condensed matter physics . Because of the explicit dependence on 
the transverse variable the dynamics of the domain wall has been reduced 
to dynamics of the collective degrees of freedom given by the Nambu-Goto 
membrane and the fields on it. In the following we shall present basic steps 
of the approach developed in papers J7|, [|. 
Step 1. The co-moving coordinate system. 

We introduce a surface S co-moving with the domain wall, that is the co- 
moving membrane. It does not have to coincide with the core, except at the 
initial instant of time. The world- volume of S, denoted by S, is parametrised 
as follows 

E9(^)( U °) = (T,r(/)). (3) 

We use the notation (w a ) a =o,i,2 = (t, cr 1 , a 2 ), where r coincides with the lab- 
oratory frame time x°, while a 1 , a 2 parametrise S at each instant of time. 
The index i = 1,2,3 refers to the spatial components of the four-vector. 
The points of the co-moving membrane S at the instant r are given by 
(F 4 )(r , cr 1 , a 2 ). The coordinate system (r, cr 1 , a 2 , £) co-moving with the do- 
main wall is defined by the formula 

x» = Y»(u a ) +£n>*(u a ), (4) 

where x M are Cartesian (laboratory frame) coordinates in Minkowski space- 
time, and (n M ) is a normalised space-like four-vector orthogonal to E in the 
covariant sense, 

n^u a )Y»{u a ) = 0, nX = 1, 

where Y% = dY^/du a . The three four-vectors Y A are tangent to S. £ is the 
transverse variable. Definition (4) implies that £ and u a are Lorentz scalars. 
In the co-moving coordinates the co-moving membrane is described by the 
condition £ = 0. For points lying on S the parameter r coincides with the 
laboratory time x°, but for £ ^ in general r is not equal to x°. 
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The extrinsic curvature coefficients K ab and induced metrics g ab on E are 
defined by the following formulas: 



K ab = n M y£, g ab = Y^ b , 



where a, b = 0, 1, 2. The covariant metric tensor in the new coordinates has 
the following form 

G a b 

1 



[Gap] 



where a, (5 — 0, 1, 2, 3; a = 3 corresponds to the £ coordinate; and 

G ab = N ac g cd N db , N ac = g ac - £K ac . 
Thus, G^ = 1, G^ a = 0. Straightforward computation gives 

v /z G = V=gh(t,if), 
where as usual g = det[g ab ], G = det[G a p], and 

fc^, « a ) = 1 - ^„ a + ^ 2 (^6 - K b a K a b ) - \eK a b K b c Kl. 

For raising and lowering the latin indices of the extrinsic curvature coefficients 
we use the induced metric tensors g ab , g ab . 

The inverse metric tensor G a/3 is given by the formula 



[G 



af3] 



G ah 
1 



where 

G ab = ( N -lyc gcd ( N -iyb_ 

N _1 is just the matrix inverse to [N ab \. It has the upper indices by definition. 

In general, the coordinates (w a ,£) are defined locally, in a vicinity of the 
world-volume E of the membrane. Roughly speaking, the allowed range of 
the £ coordinate is determined by the smaller of the two main curvature radia 
of the membrane in a local rest frame. We assume that this curvature radius 
is sufficiently large so that on the outside of the region of validity of the co- 
moving coordinates there are only exponential tails of the domain wall, that 
is that the field is exponentially close to one of the two vacuum solutions. 
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Step 2. Field equation in the co-moving coordinates. 

It is convenient to rescale the field $ and the coordinate £, 

where and s are dimensionless. We also extract from the scalar field its 
component living on the co-moving membrane and treat it separately from 
the remaining part of the scalar field. To this end we write the identity 

( p(s,u a ) = B(u a )^s) + X (s,u a ), (5) 

where 

Z?K)l0(OX) (6) 
is the component of the scalar field living on the co-moving membrane, and 

X £<p(s,u a )-B(u a )Ms) (7) 

is the remaining part. The auxiliary, fixed function ipo(s) depends on the 
variable s only. It is smooth, concentrated around s — 0, and 

y> (o) = i. (8) 

It follows that 

x (0,u a ) = 0. (9) 
The best choice for ^o( s ) is given by formula || 

1 



^o(s) 



cosh (s) 



Next we derive Euler-Lagrange equations by taking independent variations 
of B(u a ) and x- The variation 5x has to respect condition (9), hence 

5 X (0,u a ) =0. 

Because of this condition, variation of the action functional 

2 



2 f n 

S = — J dsd u y/^gh(s,u a )C 
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with respect to x gives Euler-Lagrange equation in the regions s < and 
s > 0. It has the following form 
o 1 

yab < 



M 2 



d a [^hG ab d b (B^ + X )] (10) 



+^d s [hd s (B^ + x)} + h(B^j + x)[l - m + xf] = 0. 

At s = there is no Euler-Lagrange equation corresponding to the variation 
5x- Instead, we have the condition (9). Equation (10) should be solved in the 
both regions separately, with (9) regarded as a part of boundary conditions 
for x- To complete the boundary conditions we also specify the behaviour 
of x for |^| much larger than the characteristic length l , that is for \s\ >> 1: 
we shall seek a solution such that x is exponentially close to +1 for s ^> 1 , 
while for s <C — 1 it is exponentially close to —1. 

At this stage of considerations Eq.(10) should not be extrapolated to 
s = 0. For example, the l.h.s. of it could have a 5(s)-type singularity. It 
would occur if x was smooth for s > and for s < but had a spike at 
s = 0. 

In addition to Eq.(10) we also have the Euler-Lagrange equation corre- 
sponding to variations of B(u a ). This equation has the following form 



2 r 1 



-ghG ab d b (B^ + x)\^o (11) 

-\J ds hd s ipodsWo + X) 
+ Jds hMB*Po + X) [l - (#Vo + Xf 



0. 



Here and in the following we use / ds as a shorthand for the definite integral 

Using equations (10) and (11) one can prove that x does not have the 
spike at s = 0. It follows that Eq.(10) is obeyed by x also at s = 0, and that 
now Eq.(ll) can be obtained by multiplying Eq.(10) by ^o( s ) and integrating 
over s. Therefore we may concentrate on solving Eq.(10). 
Step 3. Expansion in the width. 

First we solve Eq.(10) in the leading approximation obtained by putting 
l/M = 0. The equation is then reduced to 

^ S V°) + 0( O )[1-(0(°)) 2 ]=O, (12) 
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where 

0(o) = jB (°Vo + X (O) - 

Equation (12) does not contain derivatives with respect to time, in spite of 
the fact that it is supposed to approximate the evolution equation (10). This 
annoying fact is due to the singular character of the perturbation given in 
Eq.(10) by the terms proportional to positive powers of 1/M. We shall see 
that time evolution is obtained indirectly, from consistency conditions. 
Equation (12) has the following particular, well-known solution 

(o) =tanhs. (13) 

This solution together with conditions (8), (9) gives 

£(°)=0, x (0) =tanhs. (14) 

Notice that has the same form as the planar domain wall (2) — in this 
way we realise the idea that in the co-moving reference frame the domain 
wall does not differ much from the planar domain wall. 

The solution (13) in the co-moving coordinates does not determine the 
field in the laboratory frame because we do not know yet the position of the 
co-moving membrane with respect to the laboratory frame. Equations (10), 
(11) yield an equation for the co-moving membrane, otherwise they would 
not form the complete set of evolution equations for the field $. In fact, we 
shall see that the first order terms in Eq.(10) imply Nambu-Goto equation 
for the membrane. 

The expansion in the width has the form 

X (s, u a ) = tanh 8 + jjX (1) (s, u a ) + ^X (2) (*, u a ) + j^X (3) (s, u a ) + (15) 

B{u a ) = + ^£ (2) K) + ^ (3 V) + ( 16 ) 

where we have taken into account the zeroth order results (14). The ex- 
pansion parameter is 1/M and not 1/M 2 because 1/M in the first power 
appears in h and G ab functions after passing to the s variable. In order to 
obey the condition (9), and to ensure the proper asymptotics of x at large 
\s\ we assume that for n > 1 

x M(0X)=0, lim x W = 0. (17) 
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Inserting the perturbative Ansatz (15, 16) in Eqs.(10, 11), expanding the 
l.h.s.'s of them in powers of 1/M, and equating to zero coefficients in front of 
the powers of 1/M we obtain a sequence of linear, inhomogeneous equations 
for x {n) (s,u a ),B^(u a ) with n > 1. 

In particular, Eq.(lO) expanded in the powers of 1/M gives equations of 
the type 

L X {n) =f {n \ (18) 

where the source term is determined by the lower order terms in x an d 
B, and 

Lli9 s 2 + l-3(x (0) ) 2 . 
Explicit solution of Eq.(18) is given by the formula 

X (n) (s) = 2^i(s) f dx ^{x)f n \x) - 2^ (s) fdx ^(x)f( n \x). (19) 

J-oo JO 

This solution obeys the boundary conditions (17). 

Obviously, we assume that all proportional to positive powers of 1/M 
terms in Eqs.(10, 11) are small. For this it is not sufficient that the extrinsic 
curvatures are small, that is that IqK% <C 1. We have also to assume that 
the derivatives x^\ Bffi are n °t proportional to M. It is not the case, for 
example, if x an d B contain modes oscillating with a frequency ~ M. They 
would give positive powers of M upon differentiation with respect to u a . 
If such oscillating components were present the counting of powers of 1/M 
would no longer be so straightforward as we have assumed. This assumption 
excludes radiation modes as well as massive excitations of the domain wall. 
Therefore, the approximate solution we obtain gives what we may call the 
basic curved domain wall. To obtain more general domain wall solutions one 
would have to generalize appropriately the approximation scheme. Actually, 
the fact that such particular radiationless, unexcited curved domain wall 
exists is a prediction coming from the 1/M expansion. The expansion yields 
domain walls of concrete transverse profile — the dependence on s is explicit 
in the approximate solution we construct even at the initial instant of time. 
We may choose the initial position and velocity of points of the membrane 
but the dependence of the scalar field on the variable s at the initial time 
is given by formulas (19). This unique profile is characteristic for the basic 
curved domain wall. 
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Step 4. The consistency conditions. 

The first order terms in Eq.(lO) give the following equation 

L X m = K a a d sX {0 \ (20) 

where x^ is given by the second of formulas (14). 

The most important point in our approach is the observation that oper- 
ator L has a zero-mode, that is the normalizable solution 

Ms) = — \x~ 

cosh s 

of the homogeneous equation 

Uo = 0. (21) 

Notice that -ipo = d s x^ — this means that the zero-mode ipo is related to the 
translational invariance of Eq.(12) under s — > s + const. The presence of the 
zero-mode implies the consistency conditions. For example, let us multiply 
Eq.(20) by ip and integrate over s. It is easy to see that / tpoLx^ vanishes 
because of (21), and we obtain the following condition 

K a a Jds Ms)dsX {0) (s)=0 

which is equivalent to 

Kl = 0. (22) 

Eq.(22) coincides with the well-known Nambu-Goto equation. It determines 
the motion of the co-moving membrane, that is the functions Y l (u a ), i=l,2,3, 
once initial data are fixed. When we know these functions we can calculate 
the extrinsic curvature coefficients and the metric g a b. Review of prop- 
erties of relativistic Nambu-Goto membranes can be found in, e.g., [[|. 

Due to Nambu-Goto equation (22) the r.h.s. of Eq.(20) vanishes and the 
resulting homogeneous equation 

L X « = 

with the boundary conditions (17) has only the trivial solution 

= o. 
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Notice that vanishing x does not mean that the first order correction to 
the total field also vanishes — there is the first order contribution equal to 
B^ipo/M. It does not vanish on the co-moving membrane that is at s = 0. 

Analogous reasoning can give nontrivial consistency conditions also for 
equations (18) with n > 1. For some n, e.g. n = 2 we obtain only the trivial 
identity = 0. 

Equation (11) expanded in powers of 1/M gives equations for B^ coin- 
ciding with the consistency conditions. This follows from the fact that both 
Eq.(ll) and the consistency conditions are obtained by multiplying Eq.(10) 
by the zero-mode ipo and next integrating over s. Euler-Lagrange equation 
(11) can be regarded as generating equation for the consistency conditions. 

Equations (10), (11) in the first order do not give any restriction on the 
function . Equation for follows from the third order terms in Eq.(ll): 

^d a (^g a %B^) + (^--l)K^K b a B^ + ^-(B^) 3 = {^--l)K a b K b c K c a , 
V~9 4 35 o 

(23) 

This situation is typical for singular perturbation theories of which the 1/M 
expansion is an example — higher order equations imply restrictions (the 
consistency conditions) for the lower order contributions |10j . 
Step 5. Initial data. 

At this point we have the complete set of equations determining the evolu- 
tion of the domain wall in the 1/M expansion. Each of Eqs.(10), (11), (22) 
describes different aspect of the dynamics of the curved domain wall. Ex- 
panded in the positive powers of 1/M Eq.(10) determines dependence of \ 
on s. Because the term Bip in formula (5) has explicit dependence on s, we 
may say that Eq.(10) for \ fixes the transverse profile of the domain wall. 

Equation (11) determines the B^ n \u a ) functions, which can be regarded 
as (2+l)-dimensional scalar fields defined on X and having nontrivial non- 
linear dynamics. The extrinsic curvature K a b of S acts as an external source 
for these fields. The fields B^ can propagate along E. One may regard 
this effect as causal propagation of deformations which are introduced by the 
extrinsic curvature. 

Finally, Nambu-Goto equation (22) for the co-moving membrane deter- 
mines the evolution of the shape of the domain wall. 

Equation (22) for the co-moving membrane and equations for B^ n ' ob- 
tained from the consistency conditions are of the evolution type — we have 
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to specify initial data for them, otherwise their solutions are not unique. 
Equations (18) for the perturbative contributions x are °f different type 
- in order to ensure uniqueness of their solution it is sufficient to adopt the 
boundary conditions (17). The initial data for B(u a ) and Y l (u a ) follow from 
initial data for the original field <fi. From such data for we know the initial 
position and velocity of the core. We assume that at the initial instant r 
the co-moving membrane and the core have the same position and velocity. 
Hence, 

initial data for the membrane = initial data for the core. 

Using formula (5) one can show that then 

B^(r ,a\a 2 ) = 0, d T B^\r Ql a\a 2 ) = 0. (24) 

In order to find the domain wall solution one should first solve the col- 
lective dynamics, that is to compute evolution of the co-moving membrane 
and of the fields. The profile \ °f the domain wall is found in the next 
step from formulas (15, 19). In our perturbative scheme the profile of the 
domain wall can not be chosen arbitrarily even at the initial time — it is 
fixed uniquely once the initial data for the membrane and for the B field are 
given. Evolution of the core can be determined afterwards, from the explicit 
expression for the scalar field 0. 

* * * 



In paper || the perturbative solution has been constructed up to the 
fourth order. 

In paper |l]J second order perturbative solutions for cylindrical and spher- 
ical domain walls have been compared with numerical solutions of the Euler- 
Lagrange equation. The results of the comparison are quite encouraging. 

In paper [JT3J analogous approach has been applied to a vortex in the 
Abelian Higgs model. 

The perturbative approach to dynamics of the domain walls we have 
sketched can be generalised to models involving several fields. Also the re- 
quirement of relativistic invariance can be dropped out. 
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